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Abstract. We present dilepton spectra from nucleus-nucleus collisions at SIS energies, which 
were simulated with the GiBUU transport model in a resonance-model approach. These 
spectra are compared to the data published by the HADES collaboration. We argue that the 
interpretation of dilepton spectra at SIS energies critically depends on the couplings between 
the p meson and the baryonic resonances. 



1. Introduction 

While the vacuum properties of most hadrons are known to good accuracy nowadays, it is a 
heavily debated question how these properties change in nuclear matter. In particular, various 
theoretical predictions regarding the in-medium properties of the light vector mesons have been 
suggested. For recent reviews on in-medium physics, see [l}|3]. 

Among the expected in-medium effects, a so-called "collisional broadening" of the meson 
spectral function, due to collisions with the hadronic medium, is expected. A second class of 
predictions claims that the vector-meson masses are shifted in the medium due to the partial 
restoration of chiral symmetry |]4|. QCD sum rules can constrain these effects, but do not provide 
definitive predictions [5j. 

The more prominent hadronic decay modes of the vector mesons are unfavorable for studying 
in-medium effects, since they are affected by strong final-state interactions with the hadronic 
medium - in contrast to the rare dilepton decay modes. Since the leptons only interact 
electromagnetically, they are ideally suited to carry the in-medium information outside to the 
detector, nearly undisturbed by the hadronic medium. 

Early measurements of dilepton spectra from heavy-ion collisions in the low-energy regime 
were conducted by the DLS collaboration |6|, showing a clear excess over the expected yield. 
A similar excess was also observed in experiments at higher energies [7j|8], where it could be 
attributed to an in-medium broadening of spectral functions [9-11 . For the DLS data such 



in-medium effects never provided a convincing explanation - a problem that was soon known as 



the "DLS puzzle" 12-15 



More recently, the HADES collaboration at GSI has set up an extensive program 



for measuring dilepton spectra from p+p, p+A and A+A reactions 16-21 , in order to 
systematically check the old DLS data with improved statistics and to finally resolve the DLS 
puzzle. Up to now this endeavor has fully confirmed the validity of the DLS data and shifted 
the puzzle into the theory sector. 



2. The model 

In this paper, we apply the Gieflen Boltzmann-Uehling-Uhlenbeck transport model (GiBUU) 
[22[[23] to the nucleus-nucleus reactions studied by the HADES collaboration. The same model 
was used earlier to describe the elementary HADES and DLS data [24] . GiBUU is an hadronic 
transport model, which relies on a resonance-model description of NN and ttN collisions at low 
energies. 

In addition to the features described in [24] , the model has received a few extensions in the 
meantime: For the production of baron resonances via AW — > NR, we use angular distributions 
according to da/dt oc inspired by one-boson exchange, in order to improve upon the isotropic 
phase-space approach used earlier. At low energies, these angular distributions can become 
important due to the geometrical acceptance of the detector. In the A Dalitz decay, we now 
use the form factor from |25| , which gives only moderate deviations from the assumption of a 
constant form factor and overall seems more reasonable than the one used previously. However, 
one should keep in mind that the issue of the A form factor has still not been settled in any 
way. 

For the elementary dilepton spectra, of course only NN collisions played a role. In heavy- 
ion collisions, other types of secondary collisions can occur, such as ttN and tttt collisions. 
Details about the treatment of these processes in GiBUU can be found in [22] . For vrA^ 
Bremsstrahlung, we use the soft-photon approximation, in the same way as previously applied 
to AW Bremsstrahlung [24] . 

3. Dilepton spectra from heavy-ion collisions 

We present here preliminary results, obtained with the GiBUU transport model, of dilepton 
spectra from nucleus-nucleus collisions, namely the light C+C system at beam kinetic energies 
of 1.0 and 2.0 AGeV and the intermediate Ar+KCl system at 1.756 AGeV, all of which have 
been measured by the HADES collaboration in a fixed-target setup at the SIS18 accelerator at 
GSI. 



Table 1. Kinematic conditions of the collision systems measured by HADES and corresponding 
cuts on the single lepton momenta (all in GeV), and the maximum impact parameters in fm. 
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Table[T]gives an overview over the different reactions, including kinematic conditions, cuts and 
impact parameters used. All reactions have been simulated with impact parameters restricted 
to the range < b < 6 max Furthermore, all spectra have been filtered through the HADES 



acceptance filter [26 27 in order to account for the geometrical acceptance and resolution of 
the detector. In addition, a dilepton opening angle cut of 6 ee > 9° is applied in all cases, as 
well as the single-lepton momentum cuts listed in table [TJ matching the experimental analysis 
procedure. 

3.1. C+C 

We start by showing in ng.[l]the dilepton mass spectrum for 12 C+ 12 C collisions at beam energies 
of 1.0 and 2.0 AGeV, compared to the HADES data (which were scaled to match the simulation 
in the pion channel). As has been argued in |18|, C+C is a sufficiently light system, so that it 
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Figure 1. Dilepton mass spectrum for C+C at 1 and 2 AGeV in comparison to the data 



from 16,17 



can be regarded, in first order, as a simple superposition of elementary NN collisions, without 
significant secondary effects or modifications due to the hadronic medium. It has also been 
shown that both C+C data sets are compatible with a superposition of the NN data (when 
subtracting the r\ component, which has a different beam energy dependence) [18] . 

From this point of view, we can expect our results for the C+C collisions to reflect the quality 
of our elementary NN results. And indeed we observe an excellent description of the C+C data 
at 2 AGeV, due to the good description of pp collisions at 2.2 and 3.5 GeV 124]. The largest 
deviation in the C+C spectrum at 2 AGeV is a slight overestimation of the data in the u mass 
region, which is barely significant since the data points already have quite large error bars. Most 



other models showed a much more severe overestimation in this mass regime 16,28 



In the C+C spectrum at 1 AGeV one can observe more significant deviations, namely an 
underestimation of the data around 300 - 400 MeV. This is probably connected to the fact that 
our model's agreement to the elementary data is not as good at lower energies (in particular in 
d+p at 1.25 GeV). Improving the description of the elementary collisions would certainly also 
improve the agreement with the C+C data at 1 AGeV. 

We stress that, in contrast to other models J29 30 , the A Dalitz channel plays only a minor 
role in our cocktail, while the p production via baryon resonances is much more important. 
Fig. [T] shows separately the dominant sources of p mesons, which are given by decays of the 
Z?i3(1520) and S < 3i(1620) resonances and the process tttt — > p. We also note that the latter is 
the only production process of the p which has a physical threshold at m = 2m 7T . All other 
processes yield dilepton contributions below the 2tt threshold (but unfortunately suffer from 
poor statistics there, due to our numerical treatment). 



3.2. Ar + KCl 

The data for this reaction have been obtained by shooting 40 Ar ions on a natural KCl target. In 
the simulation we use a target of 37 Ar instead, which represents the average of the 35 C1 and 39 K 
nuclei. Fig. [2^, shows a comparison of our simulation results to the data, which clearly exhibits 
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Figure 2. Dilepton mass spectrum for Ar+KCl at 1.756 AGeV in comparison to the data 
from [19] . Left: Total spectrum; right: rj component subtracted. 



larger deviations than in the C+C case. The discrepancies mostly concern the intermediate 
mass range of 150 to 500 MeV. 

As in the C+C case at 1 AGeV, part of the missing yield might be due to disagreement with 
elementary data. However, it is unlikely that this can account for all of the missing yield, since 
our agreement with the elementary HADES and DLS data is quite good in the energy regime 
relevant for the Ar+KCl system 24 . Instead, it seems more plausible that the missing yield 



represents an actual medium effect. This is supported by the fact that the ArKCI data apparently 
shows an excess over the elementary reference cocktail [19], which is of similar magnitude as 
the discrepancy to our model. To illustrate this, fig. [2ja shows the ArKCI spectrum with the rj 
Dalitz component subtracted, compared to the NN reference spectrum (both normalized to the 
pion channel, for details see [19]). While our model underestimates the ArKCI data, it lies close 
above the elementary reference spectrum. 

The excess yield could in principle come from two distinct sources: 

(i) Dilepton emission from secondary collisions (e.g. ttN), which is in principle covered by our 
model but could be underestimated. Here also irN Bremsstrahlung could play a role, which 
is treated in our model only in soft-photon approximation (SPA). 

(ii) An in-medium modification of the spectral functions, e.g. a broadening of the p meson or 
modifications of the nucleon resonances, which are involved in the production dynamics of 
the p, most prominently the Z?i3(1520). Such spectral modifications are not included in our 
current simulation, which fully relies on vacuum spectral functions. 

In order to decide if any of these two effects can contribute to the missing yield, both should 
be studied more closely in future investigations. 



4. Conclusions 

We have shown that the GiBUU transport model provides a reasonable description of dilepton 
spectra from C+C collisions, both at 1 and 2 AGeV. At the lower beam energy, minor 
discrepancies show up, which might be connected to the underestimation of the d+p reaction 



at 1.25 GeV. At the higher beam energy of 2 AGeV, the high-mass region of the spectrum is 
well-described by a dominant contribution from p mesons produced via baryonic resonances. 
Most other models showed a significant overestimation of the data in this region. 

Regardless of the good agreement with the elementary and C+C data, our model significantly 
underestimates the Ar+KCl data in the intermediate mass region when using vacuum spectral 
functions. This excess in the data could be an indication of non-trivial effects of the hadronic 
medium and gives rise to the hope that even larger effects can be observed in the Au+Au system. 
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